Effects of Two Energy Scales in Weakly Dimerized Antiferromagnetic Quantum Spin 

Chains 
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By means of thermal expansion and specific heat measurements on the high-pressure phase of 
(VO)2P207, the effects of two energy scales of the weakly dimerized antiferromagnetic S = 1/2 
Heisenberg chain are explored. The low energy scale, given by the spin gap A, is found to manifest 
itself in a pronounced thermal expansion anomaly. A quantitative analysis, employing T-DMRG 
calculations, shows that this feature originates from changes in the magnetic entropy with respect 
to A, dS m /dA. This term, inaccessible by specific heat, is visible only in the weak-dimerization 
limit where it reflects peculiarities of the excitation spectrum and its sensitivity to variations in A. 

PACS numbers: 75.50.Ee, 75.40.Cx, 75.30Et, 75.80.+q, 71.lS.Mb, 05.10.Cc 



In recent years, many interesting and unexpected ef- 
fects have been discovered in low-dimensional antiferro- 
magnetic (afm) quantum-spin systems. Prime examples 
include the Haldane gap for spin .5=1 chains, mag- 
netization plateaux in couplcd-dimer systems, and the 
different ground states for even- and odd-leg ladder sys- 
tems, see, e.g. [l], 0, H, 0|- These phenomena reflect the 
intricate many-body quantum character of the systems 
which renders advanced mathematical methods neces- 
sary to gain a quantitative understanding of the experi- 
mental observations. In this communication, we want to 
address another non-trivial collective phenomenon, over- 
looked in the past, of the effects of two different energy 
scales caused by a weak alternation in the spin-spin inter- 
action of an S = 1/2 afm Heisenberg chain. These scales 
determine the physical properties of the interacting sys- 
tem at different temperatures. If J\ and J2 (< J\) denote 
the alternating afm coupling constants between nearest 
neighboring spins along the chain, we define Ji as the 
large scale and the spin gap A, induced by any amount 
of alternation | a < 1 with a = J2/J1, as the small 
one. The dominant feature, governing the magnetic and 
thermodynamic properties at elevated temperatures, is 
connected with J\. In the magnetic susceptibility \i f° r 
example, the maximum due to short-range order is lo- 
= 0.64Ji, practically independent of 
However as will be shown below, the small en- 
ergy scale introduces a second distinct anomaly Q , which 
can be extraordinarily large in the coefficient of thermal 
expansion. The feature observed here is no longer pro- 
portional to that in the specific heat, viz., the so-called 
Griineisen-scaling 

8] 

does not apply. 

For the study of the effects of two energy scales in the 
interesting weak-dimerization regime, the high-pressure 
variant of (VO)2P207 (HP-VOPO) appears particularly 
well suited. HP-VOPO, first synthesized by Azuma et 
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al. @, comprises one kind of S = 1/2 (V 4 +) chains with 
alternating afm Heisenberg interactions [3, [l(J, as op- 
posed to the ambient-pressure (AP) variant, where two 
slightly different chain types were identified [HI El)- No 
magnetic ordering has been found in HP-VOPO down 
to 1.5 K ~ O.OUi/fcs, the lowest temperature studied so 
far, indicating that interchain interactions are weak. 

To verify that the model of an alternating afm Heisen- 
berg chain is appropriate for HP-VOPO, we carried 
out ab initio density functional (DFT) calculations in 
the generalized gradient approximation using an LMTO 
basis. Our analysis of the DFT results in terms of 
the Nth-order muffin-tin orbital-based downfolding tech- 
nique [r3 |. shows that the dominant V-V effective hop- 
ping integrals are the intrachain hopping via the PO4 
unit, ti = 0.114 eV, followed by the nearest neighbor V-V 
intrachain hopping t% — 0.10 eV. The next non-negligible 
V-V hopping is along the &-axis, % = 0.02 eV. Since all 
other possible V-V interaction paths are more than one 
order of magnitude smaller than t\^ 1 this analysis clearly 
confirms the alternating-chain nature of HP-VOPO. 

Single crystals of HP-VOPO were prepared by trans- 
forming previously grown AP-VOPO crystals under suit- 
able pressure and temperature conditions (3 GPa/800°C) 
[lH . Structure studies reveal phase purity within the ex- 
perimental resolution of 2-3% and give parameters in ac- 
cordance with published data T4J. ESR measurements 
yield a g-factor of g = 1.977 for BLc. The magnetic sus- 
ceptibility, measured by employing a Quantum Design 
SQUID magnetometer, is very similar to the results in 
ref. [l(J. A good description of our x(T) data for BLc 
at 2-350 K is obtained by using eq. (13a) in ref. [n|, the 
above ^-factor, and the following parameters: Ji/k B = 
127K, J 2 /k B = HI K, i.e., a = 0.873 and A/k B = 31.8K, 
and a concentration of 3.9% of paramagnetic 5=1/2 mo- 
ments per mole V 4+ with an afm Weiss temperature 0\y 



2 




T[K] T[K] 



FIG. 1: Expansion coefficients of single-crystalline HP- VOPO 
along the inter-chain a- (open circles), and intra-chain c-axis 
(open squares) and for polycrystalline material a po i y (open 
triangles). Solid lines are fits discussed in the text. 



FIG. 2: Specific heat as C/T of single crystalline HP-VOPO 
(o) of mass m — 0.74 mg. Solid line corresponds to model 
calculations described in the text. 



= (2 ± 0.6) K. These parameters, being in good agree- 
ment with those given in ref . (loj , also provide a consis- 
tent description of the specific heat and thermal expan- 
sion data, see below, and indicate that impurity effects 
are of no relevance in the T range discussed here. 

The thermal expansion was measured at tempera- 
tures 1.5-200 K employing a high-resolution capacitive 
dilatometer |15|. For the specific heat measurements at 
2-30 K, a home-built AC calorimeter [16| was used. 

Figure [1] depicts the uniaxial expansivities cti(T) = 
l~ 1 dl/dT measured on an HP-VOPO single crystal along 
(i = c) and perpendicular (i — a) to the spin-chain 
direction, together with data for polycrystalline mate- 
rial a po i y . In accordance with the crystal structure, the 
ccj's are strongly anisotropic. Particularly striking is the 
anomalous in-chain c-axis expansivity a c , yielding a huge 
negative peak anomaly around 13 K followed by a change 
of sign and a maximum at 50 K. At elevated temperatures 
T > 150 K, a c becomes very small, indicating an unusu- 
ally small lattice contribution to a c . In contrast, a a and 
oipoiy show a much larger lattice expansivity with a small 
minimum at 13 K and no obvious anomaly at 50 K. 

As will be argued below, it is these two distinct fea- 
tures at 50 K and 13 K, dominating the intra-chain expan- 
sivity a c , in which the two energy scales become mani- 
fested in HP-VOPO. An indication for the low-energy 
scale, though much less strongly pronounced, is also vis- 
ible in the specific heat of HP-VOPO, shown in Fig. [2] 
The data reveal a weak variation at low temperatures, 
consistent with a spin gap, and a shoulder-like anomaly 
around 13 K. The latter feature can be discerned also in 
the C(T) data on HP-VOPO by Azuma et al. 9], where 
it has been assigned to a Schottky-type anomaly. 



For the description of the magnetic contribution C n 
to the specific heat of HP-VOPO in Fig. [21 the density 
matrix renormalization group approach for transfer ma- 
trices (T-DMRG) has been used where the free energy 
of a one-dimensional quantum system can be calculated 
by use of the Trotter-Suzuki mapping For all calcu- 
lations we have retained between 64 and 100 states and 
used the Trotter parameter e = 0.025 so that the error of 
the Trotter-Suzuki mapping is of the order 0(e 2 ) ~ 10~ 4 , 
whereas the truncation error is of the order 10" 7 . The 
phonon background C ph was approximated by a Debye 
model with a Debye temperature &d = 426 K, which pro- 
vides a reasonable description of the lattice contribution 
to the thermal expansion, see below. Instead of fitting 
the C m data with Ji and J2 as adjustable parameters, 
requiring extensive computing time, T-DMRG calcula- 
tions for Ji/ks — 127K and J^I^b = H1K were used. 
As the solid line in Fig. [5] demonstrates, this ansatz de- 
scribes the C(T) data reasonably well. Most notably, it 
proves the hump anomaly in C'(T) at low temperatures 
to be an intrinsic feature of the weakly-dimerized afm 
S = 1/2 Heisenberg chain. The departure of the solid 
line from the data at higher temperatures is assigned to 
the inadequacy of the simple Debye model in accurately 
describing the phonons in this temperature range. 

To gain more insight into the nature of the distinct low- 
temperature anomalies in HP-VOPO, a simplified model, 
denoted independent-magnon approximation (IMA) in 
the following, will be introduced, see also ref. Q. The 
significance of the IMA is verified by comparison with the 
accurate T-DMRG results. The thermodynamic calcula- 
tions performed within the IMA are based on a singlet 
ground state separated by an energy gap A from a band 
of spin-1 excitations, cf. inset to Fig.[3J To restrict the 
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FIG. 3: Magnetic specific heat divided by ksT/Ji vs. ksT/ Ji 
for the afm 5 = 1/2 alternating-exchange Heisenberg chain 
with Ji/fes = 127 K and J 2 /k B = 111K calculated by T- 
DMRG (solid line) and IMA (broken line). Inset shows one- 
magnon dispersion relations for the alternation parameter a 
= 0.873 and 0.5, calculated by using eq. (23a) in ref. [(|. 



FIG. 4: Temperature dependence of (i) dS™/d/±* (left scale) 
and (ii) dS^/dJi (right scale) for the alternating afm S — 1/2 
Heisenberg chain with Ji/fcs = 127 K and JifkB = 111 K, i.e. 
a = 0.873 and A/k B = 31.8 K. While for (i) T-DMRG results 
are compared with IMA results for different dispersion rela- 
tions given in the figure, (ii) is calculated within T-DMRG. 



maximum number of excited states to 2 , with N the 
number of S — 1/2 sites, the occupancy of a one-magnon 
state with energy E{k) and wave vector k is allowed to 
be either or 1. In zero field, where the one-magnon 
band is threefold degenerate, we obtain for the magnetic 
contribution to the free energy per spin: 

k B Td f*/ 2d ( ( E(a,k)\\ „ 

(i) 

Here d is the average spin-spin distance along the 
chains. The one-magnon dispersion E(a, k) was calcu- 
lated to high accuracy by Barnes et al. [19( . A fairly sim- 
ple parametrization, given by eq. (23a) in ref. Q , is used 
here for numerical computations of F m , the magnetic en- 
tropy S rn = ~{dF m /dT) v and C m = -T{d 2 F m / dT 2 ) v . 
Figure [3] compares the so-derived C m with results from 
T-DMRG calculations. A representation C/T is chosen 
to make the low-temperature anomaly clearer. As ex- 
pected, the IMA nicely traces the T-DMRG results at 
low temperatures, where the magnon density is small. 
In particular, it shows the peak anomaly in C m /T at 
fcsT/Ji ~ 0.1, a signature of the low-energy scale. For 
temperatures 0.06 < ksT/ ' J\ < 0.5, however, the IMA 
tends to overestimate C m by about 10% at maximum. 
The accompanied extra entropy is balanced at higher 
temperatures, so that at fc^T = 1J\ (not shown), both 
models reveal about 97 % of the full spin entropy. 

With these models at hand, a detailed description of 
the anomalous thermal expansion of HP-VOPO is pos- 
sible. If the entropy of a system can be expressed as 



a sum of separate contributions S — J2j the vol- 
ume expansion coefficient can be written as Q /3 = 
Ejft = trEj^/fllOr, with k t = -V- l {dV/dP) T 
the isothermal compressibility. This holds true also for 
anisotropic systems as long as the volume changes occur 
under hydrostatic conditions @ . For HP-VOPO, the sum 
j runs over the magnetic (j = m) and phonon (ph) contri- 
butions. The use of the Debye model and the scale invari- 
ance of S m imply that S = S ph (T/Q D ) +S m {T/J u A*), 
where A* = A/J x . With Cf = -Q D dS p h /d@ D and 
C™ = —JidS™/dJi, the phonon and magnetic molar 
specific heat contributions, respectively, one finds 

(3 = ^L(-A*T A '^+ r*C? + TV h Ct). (2) 

Vq Ol\ 

Here T A * = -dlnA*/dlnV, T Jl = -dlnJi/dlnV and 
r p/l = —dlnQo/dlnV define different Gruneisen parame- 
ters which measure the volume dependence of the respec- 
tive characteristic temperature or energy. Note that as 
a consequence of S m being a function of the two scaling 
variables T/J\ and A*, j3 (eq.[2]) contains a contribution 
oc dS™ /dA* , which has no counterpart in the magnetic 
specific heat, i.e., f3 m is no longer proportional to C m , 
viz. Griincisen-scaling does not apply. 

Equation [2] can be generalized to the uniaxial ex- 
pansivities by OLi = jfdSff/dA* + j^dS^/dJi + 
-ff h dSQ h /d<dD, i — a,b,c. The prefactors 7^ , 7/ 1 and 
jf contain the respective strain dependencies of A*, J\ 
and Qd, some combinations of elastic compliances and 
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the molar volume Vq, see [§]. Since these quantities, along 
with the Gruneisen parameters T 3 in eq. (J2]) , generally re- 
veal only a very weak variation with temperature, they 
can be treated as constants [20|. In leading order, the T 
dependence of a.; and (3 is thus given by the derivatives 
of the entropy with respect to the different energy scales. 

The two magnetic terms of interest here, computed 
numerically within the IMA (eq. and T-DMRG for 
Ji/^b — 127 K and ^/fc^ = 111K, are displayed in 
Fig. 2] as a function of temperature. For dS™/dJ\ {— 
— C™/Ji), a broad anomaly is revealed around 50 K, re- 
flecting the short-range afm correlations associated with 
Ji . In addition, dS™ jdJ\ discloses a tiny feature around 
13 K, corresponding to the low-temperature maximum in 
C m /T = -Jx/T-dSfP/dJi in Fig.El as a result of the gap 
formation. Since this anomaly is much smaller and has 
the same sign as the one at 50 K, it cannot account for the 
change in sign and the large negative peak anomaly in a c 
at 13 K. However, as shown in Fig. [4] such an anomaly 
is revealed by dS™/dA*. This feature, constituting the 
central result of this paper, reflects the sensitivity of the 
magnetic entropy to changes in the spin gap A, probed 
in a most sensitive differential way by the expansivity. 
The effect can be observed only in the weak-dimerization 
limit, where the energy scales are clearly separated. Here 
the minimum in the one-magnon dispersion is narrow so 
that the magnetic excitation spectrum is very suscepti- 
ble to small changes in A, cf. inset to Fig. [3] To simulate 
the effect of such changes on the expansivity, comparative 
calculations of 95™/9A* have been performed within the 
IMA using different dispersion relations E(a, k) delimit- 
ing the spectrum to low energies, cf. Fig.[4j While for 
E(a — 0.873, k) shown in the inset to Fig. [3J the calcula- 
tions within the IMA (broken line in Fig. [4} nicely repro- 
duce the T-DMRG results (solid line) up to the minimum 
position, a gapped cos(2ka) dispersion with the same gap 
value gives rise to an anomaly in dS™/dA* (dotted line) 
which markedly deviates from the T-DMRG data. These 
results highlight the extraordinarily high sensitivity of 
the expansivity to peculiarities of the magnetic excita- 
tion spectrum associated with the low energy scale. 

Employing the prefactors 7/ 1 , 7 A and jf and 6 d 
as adjustable parameters, the expression for the uniaxial 
expansivities can be used to fit the a a and a c data in 
Fig.Q] (solid lines). The good quality of the fit to a a for 
T > 70 K, where the phonon contribution dominates, val- 
idates the use of the Debye model with Qd = 426 K for 
satisfactorily describing the phonons in HP-VOPO. The 
lattice contribution, however, plays a minor role for a c . 
Here an excellent fit is obtained essentially by the su- 
perposition of the two magnetic contributions depicted 
in Fig. [4] properly scaled by 7/ 1 and 7 A . After having 
assured the validity of our ansatz through these fits to 
a c and a a , a modeling of the volume expansivity /3 — 
a a + on, + a c is possible by means of eq. cf. solid 
line through a po i y = (3/3 in Fig.[T] Employing a rough 



estimate for the bulk modulus of cb — k t : = (47 ± 
6) GPa from <d D = 426 K [H, the fit yields the following 
Gruneisen parameters T A = -45 ± 6, T Jl — 1.3 ± 0.2 
and T ph = 0.3 ± 0.1. Particularly striking is the large 
r A = T A * + T Jl ~ -43 ± 6, reflecting a huge pressure 
dependence of the spin gap dA/dP — A-kt-T a = -(28 ± 
8)K/GPa. A linear extrapolation to finite pressure im- 
plies that hydrostatic pressure of about (1.1 ± 0.3) GPa 
is sufficient to close the spin gap in HP-VOPO, i.e. to 
drive the system to the quantum critical uniform limit. 

In summary, the high-pressure variant of (VO)2P207 
has been used to study the effects of two distinct energy 
scales of the afm S — 1/2 weakly-dimerized Heisenberg 
chain. We show that dS™/dJ\ accounts for the signa- 
tures in the specific heat associated with short-range afm 
correlations and spin-gap formation. However, this term 
cannot explain the large low-temperature anomaly ob- 
served in the intra-chain expansivity. The latter feature 
is well described by dS™ /dA* , which probes the entropy 
contours with respect to changes in the spin gap. This 
term, inaccessible by specific heat, is observable only in 
the weak-dimerization limit and reflects the high sensitiv- 
ity of the excitation spectrum, as regards small changes 
in A. Our results, which are based on the scale invari- 
ance of the entropy S m (J 1 ,J 2 ,T) = S m (T/J 1 ,J 2 /Ji), 
can be generalized to cover the large family of magnetic 
multi-scale systems described by Heisenberg-like models. 

The authors acknowledge fruitful discussions with B. 
Liithi. The T-DMRG calculations have been supported 
by the DFG under contracts No. KL645/4-2 and GK1052. 
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